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Abstract 



A comparison of the inclusive photon energy spectrum in the radiative decay B — > X s + 7, 
measured recently by the CLEO collaboration, with the standard model is presented, using 
a S-meson wave function model and improving earlier perturbative QCD-based computations 
of the same. The dependence of the photon energy spectrum on the non-perturbative model 
parameters, pp, the 6-quark Fermi momentum in the B hadron, and m q , the spectator quark 
mass, is explicitly shown, allowing a comparison of these parameters with the ones obtained 
from the analysis of the lepton energy spectrum in semileptonic B decays. Taking into account 
present uncertainties, we estimate BR,(B — > X s + 7) = (2.55 ± 1.28) x 10 -4 in the standard 
model, assuming |Vt s |/|Vd>| = 1.0. Comparing this with the CLEO measurement BR(5 — > 
X s + 7) = (2.32 ± 0.67) x 10~ 4 implies |14|/|1/ c6 | = 1.1 ± 0.43, in agreement with the CKM 
unitarity. 



1 Supported by Schweizerischer Nationalfonds. 



1 Introduction 



Recently, the CLEO collaboration has reported the first measurement of the photon energy 
spectrum in the decay B — > X s + 7 following the measurement of the exclusive decay 
mode B — > K* + 7 reported in 1993 by the same collaboration p|. The inclusive branch- 
ing ratio and the photon energy spectrum allow a less model-dependent comparison with 
the underlying theory, more specifically the standard model (SM), as compared to the ex- 
clusive decay modes which require additionally decay form factors. The CLEO data have 
been compared in [|IJ with the SM-based theoretical computations presented in - ||, allow- 
ing to draw the conclusion that agreement between theoretical predictions and experiment is 
good, given that large uncertainties exist in both. In particular, the measured branching ratio 
BR(5 -»• X s + 7) = (2.32 ± 0.67) x KT 4 is in agreement with the SM-based estimates in 
[§], |5|, as well as with the ones in 0, [/J. 

In this letter, we would like to report on an improved calculation of the photon energy 
spectrum compared to what we have presented earlier and which has been used in the CLEO 
analysis |l| . The main theoretical difference lies in the inclusion of the complete operator basis 
Oi, ...,08 for the effective Hamiltonian, defined below, in the computation of the partonic pro- 
cesses b — > s + 7 and b — > s + g + 7, and in using the complete leading- logarithmic computations 
of the anomalous dimension matrix presented in |§. In contrast, our previous calculations were 
done in the truncated approximation, where we had dropped the effects of the four-Fermi oper- 
ators 3 , ...,06 an d the chromomagnetic operator 8 in the computation of the contributions 
from b — > s + g + '-f. In addition, use was made of the anomalous dimension matrix derived in B. 



The other ingredient of our calculation, namely a specific I?-meson wave function model [|1(], [11 
to incorporate the non-perturbative effects on spectra, remains unaltered. However, since data 
are now available, we fit the normalized photon energy spectrum with the improved theoretical 
framework to determine from the shape the non-perturbative parameters of the wave function 
model being used, namely the Fermi motion parameter pp and the spectator quark mass m q . 
These in turn determine within the model the 6-quark mass. There is considerable theoretical 
interest in these parameters, in particular pp, which is a good measure of the kinetic energy of 
the b quark in the B meson, and the 6-quark mass. Since a similar framework has also been 
used, in conjunction with the perturbative QCD-improved parton model, in the analysis of the 
lepton energy spectrum in inclusive 5-meson decays 0, it is interesting to compare the model 
parameters obtained from the two measurements. While, admittedly, present errors are large 
preventing us from drawing sharp conclusions, some valuable insight on the shape parameters 
and normalization can already be obtained and we quantify this information. 

This letter is organized as follows: In section 2, we briefly summarize the effective Hamil- 
tonian for the decay B — > X s + 7 and present the Wilson coefficients numerically. Section 3 
contains an anatomy of the partonic processes 6^5 + 7 and b — > s + g + 7, where the essential 
steps in the derivation of the matrix elements are given. The photon energy spectrum at the 
partonic level is derived in section 4, including a discussion of the Sudakov behaviour in the 



end-point region. Section 5 summarizes the I?-meson wave function model [10], |1TJ. Numerical 
results for the branching ratio BK(B — > X s + 7) in the SM, the ratio of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements | Vd,|, and fits of the CLEO photon energy spectrum from 
B —>■ X s + 7 with our theoretical estimates, yielding a ±lcr-contour in the parameter space 
(pF,m q ) are given in section 6. 
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2 Effective Hamiltonian for the decay B — > X s + 7 



The framework we use here is that of an effective theory with five quarks, obtained by integrating 
out the heavier degrees of freedom, which in the standard model are the top quark and the 
PU-boson. A complete set of dimension-6 operators relevant for the processes 5 ^ s + 7 and 
6^s + 7 + gis contained in the effective Hamiltonian 



H eff (b -> ary) 



AC 8 

v2 j= i 



(1) 



where Gf is the Fermi constant coupling constant and Cj(/i) are the Wilson coefficients eval- 
uated at the scale /z, and X t = V t bV t * s with Vij being the CKM matrix elements. The operators 
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(2) 



where e and g s are the electromagnetic and the strong coupling constants, respectively. In the 
magnetic moment type operators O7 and Os, F^ u and G^ v denote the electromagnetic and the 
gluonic field strength tensors, respectively. L = (1 — 7s)/2 and R = (1 + 7s)/2 denote the 
left and right-handed projection operators. We note here that the explicit mass factors in O7 
and 0§ are the running quark masses. QCD corrections to the decay rate for b — > S7 bring in 
large logarithms of the form a™ {my/) log m (m^/M), where M = m t or my/ and m < n (with 
n = 0, 1,2, ...). To get a reasonable result, one should resum at least the leading logarithmic 
contribution (i.e. m — n) to all orders. Using the renormalization group equation the Wilson 
coefficient can be calculated at the scale \i ~ m& which is the relevant scale for B decays. At 
this scale the large logarithms are contained in the Wilson coefficients. As we are working in 
this paper to leading logarithmic precision, it is sufficient to know the leading order anomalous 
dimension matrix and the matching Ci(fi = my/) to lowest order (i.e., without QCD corrections) 
The 8x8 matrix is given in || and the Wilson coefficients are explicitly listed in || [13|| . 
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The numerical values of the Wilson coefficients that we use in our calculations are given in table 
1. For subsequent discussion we define two effective Wilson coefficients C7 (//) and Cg (//) 
below and give their numerical values in table 1: 



c e J f 



= C 7 + QdC 5 

= c 8 + c 5 . 



3Qd Cq 



(3) 



2 The results given here for the entries concerning O7 and Os correspond to the naive dimensional regular- 
ization scheme (NDR), which we use in the calculation of all the matrix elements. 
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Ci(/i) 


/X = mjy 


fj, = 10.0 GeV 


ji = 5.0 GeV 


fj, = 2.5 GeV 


Ci 


0.0 


-0.158 


-0.235 


-0.338 


c 2 


1.0 


1.063 


1.100 


1.156 


c 3 


0.0 


0.007 


0.011 


0.016 


c 4 


0.0 


-0.017 


-0.024 


-0.034 




0.0 


0.005 


0.007 


0.009 




0.0 


-0.019 


-0.029 


-0.044 


c 7 


-0.193 


-0.290 


-0.333 


-0.388 




-0.096 


-0.138 


-0.153 


-0.171 




-0.193 


-0.273 


-0.306 


-0.347 


r< e ff 


-0.096 


-0.132 


-0.146 


-0.162 



Table 1: Wilson coefficients Cj(/i) at the scale /i = m w = 80.33 GeV ("matching conditions") 
and at three other scales, \i = 10.0 GeV, fi = 5.0 GeV and /i = 2.5 GeV, evaluated with two- 
loop /3-function and the leading-order anomalous- dimension matrix. The entries correspond 
to the top quark mass rnt{rr?° U l ) = 170 GeV (equivalents, m p t ole = 180 GeV) and the QCD 
parameter with 5 flavours A 5 = 195 MeV (equivalently, a s (m 2 z ) = 0.117), both in the MS 
scheme. 



3 An anatomy of the processes b — > 57 and b —> s^g 

We summarize the principal points of the derivation of the photon energy spectrum in the decay 
b — > s r y(+g) in this letter and refer to for technical details. As it is the process b — > s^g 
which leads to a nontrivial photon energy spectrum at the partonic level, we strive at taking it 
fully into account. The amplitude for this decay suffers from infrared singularities as the gluon 
or the photon energy goes to zero. The topology of the decay b — > s'yg resembles that of the 
two-body decays b —> 57 and b — > sg in the limit as E g — > or — > 0, respectively. The 
singular configurations in b — > sjg are cancelled in a distribution sense in the photon energy 
spectrum if one also takes into account virtual corrections to the two-body process b — > sg and 
b — > 57, order by order in perturbation theory. We will take into account only those virtual 
correction diagrams which are needed to cancel the infrared singularities from b — ► 57^. 

We first discuss the result for b — > sj. It turns out that the effects of the four- Fermi 
operators to b —>■ 57 can be absorbed into a redefinition of the coefficient C-j — > C7 , with C7 
defined in eq. ([3]). This not only holds for the b —>■ sj amplitude without virtual corrections 
but also for those infrared- sensitive virtual corrections which we have mentioned above. One 
therefore has to calculate the matrix element of C; O7 for b — > 57 including virtual corrections. 
The result of this calculation, which was derived in d = 4 — 2e dimensions, can be expressed as 

yvirt -pvirt , -pvirt (a\ 

1 7 ~ L 7,sing ' 1 7 Jin ) V 1 ) 

where we have split this quantity into an infrared- finite and an infrared- singular part: 
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(log 2 r — 4 log r + 4 logr log(l 



3 log r + 8 log(l — r) — 2 log — g 



The quantity r is denned as r = (m s /m b ) 2 . The infrared singularity (|5|) will be cancelled when 
taking into account the gluon bremsstrahlung diagrams involving the operator O-j. 

As we already noted, the process b — > s^g has also an infrared singularity as the photon 
becomes soft. These singularities are cancelled analogously by virtual photon corrections to 



the matrix element for b —>■ sg, i.e. in Cg-^Og. The result is 
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(9) 
(10) 



where r is given in eq. (0) and = — 1/3. A remark concerning the quark masses is in order 
here. When calculating the matrix elements we have used the on-shell subtraction prescription 
for the quark masses. Due to the explicit factors of the running quark masses in the operators 
O7 and 08, the m\ factor contained in Y™ rt and Y™ rt given above should be replaced by the 
following product 



ml 



ml 



% P oie m b(v) , (11) 

where m b)PO i e and mj(ju) denote the pole mass and the running mass of the b quark, respectively. 
In actual practice, we identify all the masses m;, in the various intermediate expressions with 
mib, P oie and multiply at the end the so-derived decay width T(B — > X s + 7) with a correction 
factor R: 

R= (m b (fi)/m btPO i e ) 2 , (12) 

as also advocated in 0. 

We now take up the matrix elements for the process b — > s^g. As the explicit expressions 
are too long to be presented here, we only point out the basic structure and give the complete 
formulae elsewhere |14j| . We first concentrate on the contributions of the four- Fermi operators 
Oi, ...,Oq. It turns out that the diagrams which do not involve both gauge particle radiation 
from an internal quark are either zero or can be absorbed into a redefinition of the coefficients 
(CV, Cg) — ► (Cy , C| ). The remaining case, in which both gauge particles are emitted from 
the internal fermion line, is discussed now. There are two such diagrams associated with each 
four-Fermi operator, whose sum is ultraviolet (and infrared) finite. We denote these matrix 
elements by Mi 

4-G 

Mi = —j=-(s^g\CiOi\b) ,(z = l,2,...,6) . (13) 

Analogously, the matrix elements of C7 O7 and C% Og are denoted by M 7 and M 8 , respec- 
tively. Adding all these contributions, the complete matrix element for b s^g is denoted by 



(7) 
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M 

8 

M = Y,M l . (14) 

i 

When squaring M and summing over the polarizations and spins of the particles, it turns out 
that only |M 7 || and |M 8 || are infrared- singular. The other squared amplitudes and all the 
interference terms are infrared-safe. We therefore make the decomposition 

\M\l = F+\M 7 \l + \M 8 \l , (15) 

where F denotes the infrared-safe contributions. 



4 The photon energy spectrum at the partonic level 

In the following, it is useful to define the dimensionless photon energy x through the relation 

1 2m b 

x then varies in the interval [0, 1]. The contribution of the F-term in eq. flT5D to the spectrum 
dTp/dx is completely finite and is worked out numerically. As the contribution dT h 7 rems /dx 
associated with |My|| is singular, we work it out analytically in d = 4 — 2e dimensional phase 
space. Making use of the abbreviations r = (m s /mj,) 2 and £ = (1 —r)x the quantity dT\ rems / dx 
can be written as 

d^r ems rn w 4£ 3T 2£ r M | r I x(2x 2 - 5x - - r) 



dx £V (1 - s) 1+2£ 1 1 - £ 

(1 — r)x(l — x)(2x — 1) , . , , . .1 , . 

+ ~ (1 _ a L 2(l + i)bg(l-0-U^) . (16) 



with 
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4r 4 

4 = -4- — --(l + r)\og(l-0i 

h ^ ^log(l-0-^log(l-0-^ + 8(l + ^-^log 2 (l- . 

(17) 

Here the symbol Li stands for the Spence function. The photon energy spectrum away from 
the endpoint x = 1 can be obtained by taking the limit e — > in eqs. ([16]) and (|TT|). However, 
as the experimentally interesting region is just near this endpoint one has to add the virtual 
QCD corrections to the two-body process b — > sj, as already mentioned. Before we discuss this 
further, we give the result for dY h % ems /dx. 

As dT^~ ems /dx has a non-integrable singularity at x = (i.e., E 7 — 0), we should work in 
ci-dimensions as well. However, as the photon energy spectrum is of no experimental interest 



5 



at very small energies, i.e. — > 0, we remove the infrared regularization immediately. As 
for the branching ratio, we note that the dimensionally regularized total integral r| rems can be 
obtained from the analogous expression for Y\ rems by doing obvious replacements. The result 
for dTg rems /dx reads 

j-pbrerns \r~i \ r\ /^ e ff\2 

— = mt ; a em a s 1 + r 1 — r 1 — x) 

dx b 96tt 5 V A ; V ; 



x [ 2[(l-rMx-2) + 2(l + r)] 1 - £ | 

{ x(l — x)(l — r) r 

(l-r)(l-x)(l-2x) (l-r)(2x 2 -a; + l) 8' 
(1-0 2 W 



(18) 



Returning to the discussion of the infrared singularity of the quantity dT 7 rems /dx in eq. fll6D, 
which occurs for x — *■ 1, i.e., in the experimentally interesting region, we recall that this singu- 
larity is cancelled in a distribution sense if one takes into account the virtual QCD corrections 
to the tree level matrix element of O-j for the two-body process b — > 57, given in eqs. (|5|) and 
(H). Technically, it is useful to define the integrated quantity 



r 7 (s ) 



J'nbrems 



dx , (19) 



in which the singularities cancel manifestly. The expressions for dT\ rems / dx and Tj trt are given 
in eqs. fljjD and (f|), respectively. The explicit expression for r 7 (s ) reads 

r 7 ( So ) = v (1 + -i n) e(i - 5 ) , (20) 

OTT 

where 6 is the step function and V is 

y = mg |GF 3 A ff /|2 a em (l + r)(l-r) 3 . (21) 

The expressions for Q is identical with the one that we had derived and presented earlier [3j]. The 
difference lies in the normalization factor V given above. We have now included the complete 
operator basis, Oi, 0$, as opposed to our earlier calculations || [§] where use was made of a 
truncated basis neglecting the four-Fermi operators O3, ■■■,0^. The difference reflects itself in 
the coefficient C7 used in P, f|, which is now replaced by C% . The effective coefficient C7 
calculated with the anomalous dimension matrix given above is typically 10% smaller than C-j 
used by us earlier. From the definition of the quantity r 7 (so) it is clear that the photon energy 
spectrum is reconstructed by differentiation, i.e., 

dr^x) _ dY 7 (s ) 
dx {X) ~ ds lso=x ■ [ZZ) 

The end-point spectra, however, show sensitivity to the leftover effects of the (cancelled) in- 
frared singularity, with the photon-energy distribution rising very steeply near the end-point, 
x ~ 1. Furthermore, there is still a ^-function contribution sitting at x = 1, stemming from 
differentiation of the G function in eq. (p0|). To remedy this, one often resorts to (an all order) 
resummation of the leading (infrared) logarithms. This resummation is done at the level of 
the quantity r 7 (so), i.e., before one derives the spectrum by the differentiation just described. 
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Although we are using for Q the expression for a non-zero strange quark mass (i.e. r 7^ 0) in 
the numerics [§[], it is nevertheless instructive to look at the limit r — > 0. Splitting Q into two 
parts Q = + Q 2 and using the notation si = (1 — s ) we get 



fii = -2 log si - 71ogsi , 

u 2 4tt 2 2 
n 2 = -21og-^2 - — -5 + 10si + s 2 --^ + si(si-4)logsi. (23) 



Before giving the exponentiated analogue of eq. (|20|) , we point out that the lepton energy 
spectrum in the inclusive semileptonic decays b — > u(+g)£u£, in the limit of neglecting the final- 
state quark mass, has a similar structure |L5], [16). The leading term, i.e. the Sudakov double-log 
flTTf in eq. (|2"3]), also enters in the lepton energy spectrum in the mentioned semileptonic decays 
with the same coefficient. However, as shown here and earlier [|3], || (and also discussed in 
|i~8|), the coefficients of the non-leading logarithmic and power terms in the inclusive decay 
B —>■ X s + 7 and the semileptonic decay B — > X u £v£ decay are different []. Therefore, the 
photon energy spectrum in B — ► X s + 7 at the parton level can not be gotten from the lepton 
energy spectrum in the decays B — > X u lvi. For a recent discussion of the Sudakov- improved 
photon energy spectrum in B — > X s + 7, see |l9| . Finally, the exponentiated analogue of r7(s ) 
in eq. (BCT) reads 



r-»(a ) = V (l + ^ fi 2 ) exp 0(1 - So ) . (24) 

The explicit expressions for fli and fl 2 for r 7^ are given in [|J . The expression for the photon 
energy spectrum can again be obtained by differentiation: 



dT 7 p (x) _ dT 7 p (s ) 
dx ds s °~ x 



(25) 



As exponentiation is required only near the end-point x — > 1, we use the exponentiated form in 
the region x > x cr u only. For numerical calculations we take x cr i t = 0.85, where the transition 
from one form to the other works smoothly. 

To summarize, the final answer for the photon energy spectrum for the process b — > s r y(+g), 
where the b quark decays at rest, is given by 

dT dT F nt . dT b 7 rems n . . dT e 7 xp dT% ems . . 

-r = -r- + ®(x C rit - x) — + 6 (x - x crit ) + —5 — , 26 

dx dx dx dx dx 



where dTp/dx is too long to be given here and is relegated to |TJ|]; dT? ems /dx, dT 7 xp /dx and 



dTg ms /dx are given in eqs. (|i~6"l), ( p5|) and (|I8|), respectively. A remark is in order here: As it 
stands, eq. ( PE| ) has still a singularity at the experimentally uninteresting point £7 7 = 0, which 
is cancelled in a distribution sense when taking into account the virtual photon correction (see 
eq. (HD) of the operator Og. In principle, one could do a resummation of the soft photons, in 
analogy with the treatment of the resummation of soft gluons just discussed above. As it turns 
out numerically, the contribution of photon energies below 50 MeV to the branching ratio is 
less than 1%, we just make a cut there and do not resort to any resummation in the E 1 — > 
limit . 



3 The discrepancy in the analytic results for the non-leading terms in fig] and |16| for the semileptonic decay 
seems to have been settled in favour of lllql . 



7 



5 B-meson Wave function Effects in B — > X s + 7 

In order to implement the 5-meson bound state effects on the photon energy spectrum, we 



continue to use the wave-function model UTOl llll that we have used in our earlier work 



We recall that this model has also been used in calculating the lepton energy spectra in inclusive 
D- and i?-semileptonic decays. Assuming that the 5-meson wave function effects are universal, 
the parameters of this model can be fixed from the semileptonic analysis and one can make 
a parameter-free comparison of the photon energy spectrum in B — > X s + 7 with data. This 
is what has been done in the experimental analysis of the data However, in general, the 
non-perturbative effects are likely to depend on the final-state quark mass, which, for example, 
is the case in the decays B — > X s + 7 and B — > X c iv?. It is, therefore, an interesting question 
to ask if the non-perturbative model being used by us consistently describes the lepton and 
photon energy spectra in the inclusive decays B — ► Xlp^ and B —>■ X s 7, respectively. To answer 
this question quantitatively, one has to implement the perturbative QCD effects, as discussed 
in the previous section, incorporate the model, and confront the resulting framework with data 
to determine the model parameters with a well-defined statistical significance. This is what we 
shall undertake in this paper. We note that this point has also been investigated recently with 
a somewhat different non-perturbative (model) and perturbative(QCD) input in |HJ. Due to 
detailed differences in the underlying theoretical frameworks and in the analysis of data, no 
direct quantitative comparison with this work is attempted here. 

In the present model, which we refer to as the Fermi motion model, the 5-meson consists 
of a 6-quark and a spectator q and the four-momenta of the constituents are required to add 
up to the four-momentum of the i?-meson. In the rest frame of the 5-meson the 6-quark and 
the spectator fly back-to-back with three momenta pb = —p q = p. Energy conservation then 
implies the equation 



m B = ymf + fp + ^m 2 q + fp , 

which can only hold for all values of \p\, if at least one of the masses becomes momentum 
dependent. In this model the spectator quark m q is chosen to be a fixed, i.e. momentum- 
independent, parameter and the 6-quark mass becomes momentum dependent. This momentum 
dependent 6-quark mass is denoted by W in the following and is given by 



W 2 (p) = M B 2 + m q 2 -2M B ^/p 2 + m q 2 . (27) 

The 6-quark, whose decays determine the dynamics, is given a non-zero momentum having a 
Gaussian distribution, with the width determined by the parameter pp: 

0(p) = ^i^exp(^) ; p=\p\ (28) 

y^PF PF 2 

We note that this wave function is normalized such that 

dpp 2 4>{p) = 1 



00 







The photon energy spectrum from the decay of the 5-meson at rest is then given by 

(IT fPmax dY h 



dEl L ^m^V^) , (29) 

where p m ax is the maximally allowed value of p and 4tt- is the photon energy spectrum from 
the decay of the 6-quark in flight, having a mass W(p) and momentum p. This can be obtained 
by Lorentz boosting the 6-quark decay spectrum at rest, which has been obtained in eq. (p6|). 
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6 Estimates of BR(5 -> X s + 7) in the SM and the Pa- 
rameters (pi? : m q ] 



It is theoretically preferable to calculate the branching ratio for the inclusive decay B —>■ X s + , ~f 
in terms of the semileptonic decay branching ratio 

BR(B -> X sl ) = [ r(£? ^ 7 + * s) p BR(5 - X^) , (30) 

J- si 

where, in the approximation of including the leading-order QCD correction, T s i is given by the 
expression 

r " = G loo C f m 5 b g(m c /m b ) (1 - 2/3^/(m e /m»)) . (31) 
The phase space function g(z) is defined as 

g(z) = l-8z 2 + 8z 6 - z 8 - 2Az 4 ln(z) , (32) 
and the function f(z) can be seen, for example, in ref. ||10[| . An approximate analytic form for 



f(z) has been obtained in [15 



/W = (vr 2 -y)(l-^) 2 + ^ • (33) 



To get the branching ratio in eq. ( |30"D we can proceed in two ways. One can either take 
the partonic (purely perturbative) expressions for both T(B — > X s 7) and or one can first 
implement the wave function effects and then integrate the spectra. In the latter case, the 



dominant wave function effects in the quantity T s i are included if one identifies m& in eq. (31) 
with the effective value W e ff of the b quark mass, which is the value of the floating 6-quark 
mass averaged over the Gaussian distribution: 

W* /f = J dpp 2 ^p)W 5 (p) , (34) 

where W(p) is given in eq. (p7|) . We remark that these procedures yield an almost identical 
branching ratio (within 1%), which shows that the Gaussian-distributed Fermi motion model 
[|i~0| , |i~T|1 is in agreement with the result that power corrections to the inclusive decay widths 
T(B — > X s + '~f) and T(B — > Xivt), calculated in the context of the heavy quark effective theory 
[ p0|| , cancel out in their ratio. 

We now estimate BK(B —>■ X s + 7) in the standard model and theoretical uncertainties 
on this quantity. The parameters that we have used in estimating the inclusive rates for 
BR(5 —>■ X s + 7) are summarized in table |2|. The largest theoretical uncertainty stems from 
the scale dependence of the Wilson coefficients as was already stressed in [^TJ. The extent of 
this variation is somewhat correlated with the value of the QCD scale parameter A5 and the 
top quark mass. As given explicitly in the preceding section, the decay rate for B — > X s + 7 
depends on seven of the eight Wilson coefficients given earlier, once one takes into account the 
bremsstrahlung corrections and is not factored in terms of a single (effective) coefficient, that 
one encounters for the two-body decays b — > s + 7 0. To get some insight in the errors we 
enumerate here the values of the two dominant effective coefficients, Cy and , as one 
varies /i, A 5 and m t in the range given in table 2: 

Cj ff = -0.306 ±0.050 

C e J f = -0.146 ±0.020 (35) 
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Parameter 


Range 


rn t (GeV) 
\i (GeV) 
A 5 (GeV) 
B{B -> X£z/ £ ) 
m c /m b 
m w (GeV) 

a QED 


170 ±11 

r n+5.0 
5.0 25 

0.195lgr 
(10.4 ±0.4)% 
0.29 ±0.02 
80.33 
130.0 



Table 2: Values of the parameters used in estimating the branching ratio BK{B — > X s + 7) in 
the standard model. The range of A5 is taken from the present world average (corresponding 
to a s (m%) = 0.117 ± 0.005, using the two-loop /3-function p2| ) and the semileptonic branching 
ratio from 



The present theoretical uncertainties on these coefficients represent the dominant contribution 
to the theoretical error on BR(£> — > X s + 7), contributing about ±35%. The second source of 
scale-dependence is due to the appearance of the running quark masses in the operators O7 and 
0$. As discussed in section 3, this brings into fore the extra (scale-dependent) multiplicative 
factor R = (ra&(//) / m^poie) 2 for the branching ratio BR(B — > X s + 7). The next largest error 
arises from the parameters which are extrinsic to the decay B — ► X s + 7 but have crept in 
due to our procedure of normalizing the branching ratio BR(i? — > X s + 7) in terms of the 
semileptonic branching ratio. To estimate this extrinsic error, we note that the ratio m c /mb 
can be written as •m c jm b = 1 — (mj, — m c )/m&. The mass difference m& — m c is known to a 
very high accuracy through the 1/toq expansion or from the semileptonic b — > c spectrum; 
to a rather high accuracy its value is — m c = 3.40 GeV p5| . The b-quark mass m& is, 
however, not known so precisely. Using for the b quark pole mass m b)PO i e = 4.8 ± 0.15 GeV, 
one gets m c jm b = 0.29 ± 0.02, which is consistent with the determination of the same from 
the lepton energy spectrum m c /rrib = 0.316 ± 0.013 but less precise. This introduces 
considerable uncertainty in the theoretical estimates of the branching ratio BR(I? — > X s + 7), 
which can be judged from the numerical value g(z = 0.29 ± 0.02) = 0.542 ± 0.045, leading to 
±8.1% error on the branching ratio BK(B —>■ X s + 7). We recall here that f(z) is a slowly 
varying function of z, and for the quark mass ratio relevant for the decay b —>■ c+ £ namely 
z = 0.29 ± 0.02, it has the value f(z) = 2.57 =F 0.06. Taking into account the experimental 
error of ±4.1% on BR(i? — > Xivg) p3]| , and adding these errors linearly, one gets an error of 



±12% on BR(£> — > X s ±7) from the semileptonic decays. The procedure that we have adopted 
in estimating the theoretical uncertainties on the inclusive branching ratio BR(Z? — > X s ± 7) 
is as follows: We propagate the errors due to the scale /x, the QCD scale parameter, A 5 , 
and the top quark mass in our calculations. As remarked already, this constitutes the largest 
error. The extrinsic errors (from m c /m;, and the semileptonic branching ratio), being obviously 
uncorrelated, are then added in quadrature in quoting the branching ratio. 



We now proceed to discuss our results. Assuming |Vts|/|Kb = 1 HU, we Pfo^ m Fig- 1 the 
branching ratio BR(5 — > Xg'j) as a function of the top quark mass m t . For all three solid 
curves, representing the SM-branching ratio, we have used z = 0.29. The top solid curve is 
drawn for /i = 2.5 GeV and A 5 = 0.260 GeV. The bottom solid curve is for fi = 10 GeV and 
A 5 = 0.145 GeV, and the middle solid curve corresponds to the central values of the input 
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Figure 1: BR(i? — > X s 7) as a function of m t . The three solid lines correspond to the variation 
of the parameters (/x, A 5 ) as descibed in the text. The experimental (±lcr)-bounds from CLEO 
|I| are shown by the dashed lines. 

parameters in table 2. Using m t = (170 ± 11) GeV, and adding the extrinsic error, we get 

BR(B -f X s + 7) = (2.55 ± 1.28) x 10~ 4 , (36) 

to be compared with the CLEO measurement BR(B -> X s + 7) = (2.32 ± 0.67) x 10" 4 . The 
(±l<j)-upper and -lower bound from the CLEO measurement are shown in Fig. 1 by dashed 
lines. We see that the agreement between SM and experiment is good. The theoretical errors 
estimated by us are, however, larger, than for example in |J, for reasons that we have explained 
above. 

In Fig. 2 we show the branching ratio BK(B — > X s + 7) as a function of the CKM matrix 
element ratio squared | V ts \/\ V c b\ 2 , varying m t , \i and A 5 in the range specified in table 2. Using 




1 2 

(IVJ/IVJ) 2 



Figure 2: BR(£> — ► X s j) as a function of (| Vj s |/| V^,|) 2 . The solid lines correspond to the 
variation of the parametrs (fi,A^,,m t ) in the limits specified in table 2. The experimental 
(±lcr)-bounds from CLEO |l|] are shown by the dashed lines. 



the (±lcr)-experimental bounds on the branching ratio (dashed lines) we infer 

\V ts \/\V cb \ = 1.10 ±0.43, 



(37) 
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-Ely-interval 


number of events 


1.95 - 2.20 GeV 
2.20 - 2.45 GeV 
2.45 - 2.70 GeV 
2.70 - 2.95 GeV 


229 ± 256 
484 ± 163 
381 ± 105 
12 ±59 



Table 3: Photon yield in the laboratory frame from the decay B — > X s + 7, obtained from 
the measurement of the photon energy spectrum by the CLEO collaboration based on a 
sample of 2.152 million BB events. The data shown have been corrected due to the detector 



acceptance j2"T | 



which is consistent with the indirect constraints from the CKM unitarity |[22| . 

Now we discuss the photon energy spectrum and fit of the parameters pf and m q by using 
the CLEO data which has been corrected for the detector effects. The photon yield from the 
decay B — > X s + 7 is given in table 3 in photon energy bins having a width of 250 MeV starting 
with E 1 = 1.95 GeV. We note that these entries are based on the weighted average of the 
two different methods (event shape and B reconstruction) used by the CLEO collaboration 
in the analysis of their data P7| . The photon energy has been measured in the laboratory 
frame (i.e. in the rest frame of T(4S')) and the numbers in table 3 are presented in this frame. 
This implies that the B mesons from the T(4S") decay have a momentum of ~ 350 MeV, and 
in doing the analysis we have boosted the theoretical rest-frame spectra accordingly. As the 
theoretical uncertainties are mainly in the normalization of the spectrum, we normalized both 
the theoretical predictions (parametrized in terms of pf and m q ) and the experimental data to 
unit area in the interval between 1.95 GeV and 2.95 GeV. We then performed a \ 2 analysis. The 
experimental errors are still large and the fits result in relatively small \ 2 values; the minimum, 
Xmin = 0.038, is obtained for p F = 450 MeV and m q = 0, which corresponds to the 6-quark pole 
mass m^poie = 4.77 GeV, in good agreement with theoretical estimates of the same. A contour 
plot in the (m q ,pF) plane, obtained by varying the \ 2 by one unit from xLw which corresponds 
to ~ 39% confidence level , is shown in Fig. 3. As expected, due to the large errors in data 



1000 




300 400 
Fig. 3 



Figure 3: Contourplot in the (m q ,pF) parameter space obtained from \ 2 — Xmin + IThe 
minimum \ 2 is for the values (0 MeV, 450 MeV). 
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the parameter space cannot be restricted very much at present. Nevertheless, the results for 
Pf are compatible with the value Pf = 270 ± 40 MeV obtained from the 5-semileptonic decay 
analysis [ffj. In Fig. 4 we have plotted the photon energy spectrum normalized to unit area 
in the interval between 1.95 GeV and 2.95 GeV for the parameters which correspond to the 
minimum x 2 (solid curve) and for another set of parameters that lies near the x 2 -boundary in 
the contour plot (dashed curve). Data from CLEO Jl| are also shown. 




E 7 [GeV] Fig . 4 

Figure 4: Comparison of the normalized photon energy distribution using the corrected CLEO 
data |1J and our theoretical distributions, both normalized to unit area in the photon energy 
interval between 1.95 GeV and 2.95 GeV. The solid curve corresponds to the values with the 
minimum x 2 , (m q , pi?) = (0,450 MeV), and the dashed curve to the values (m g ,pi?) = (300 MeV, 
310 MeV). 

In Fig. 5 we show the comparison of the differential branching ratio from CLEO with our 
calculations. The theoretical curves correspond to the central values of the input parameters 
in table 2. The agreement between experiment and SM is good. 




2 2.2 2.4 2.6 2.8 3 



E,[GeV] Fig . 5 

Figure 5: The same as in Fig 4, but with absolute (i.e not normalized) differential branching 
ratio. 

In summary, we have presented an improved theoretical calculation of the branching ratio 
BR(5 — > X s + '-f) and the photon energy spectrum in B — > X s + / y, using perturbative QCD and 
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a Gaussian Fermi motion model. We estimate BR(B -> X s +^) = (2.55±1.28) x 10~ 4 in the SM, 
in agreement with the corresponding branching ratio BK(B — > X s + 7) = (2.32 ± 0.67) x 10 -4 
reported by the CLEO collaboration. Our model calculations provide a good account of the 
measured photon energy distribution, and the best-fit parameters correspond to pp = 450 MeV 
and m^poie = 4.77 GeV. The errors on these parameters are still large but within errors both 
p F and m bjPO i e determined from the radiative and semileptonic B decays are compatible with 
each other. Precise comparison requires improved measurements and theory, which we hope 
are forthcoming. 
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